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A hidden order that emerges in the frustrated pyrochlore Tb 2 +j:Ti 2 -j: 07 + 3 , with Tc = 0.53 K is studied using spe¬ 
cific heat, magnetization, and neutron scattering experiments on a high-quality single crystal. Semi-quantitative 
analyses based on a pseudospin-1/2 Hamiltonian for ionic non-Kramers magnetic doublets demonstrate that it is 
an ordered state of electric quadrupole moments. The elusive spin liquid state of the nominal Tb 2 Ti 207 is most 
likely a U(l) quantum spin-liquid state. 

PACS numbers: 75.40.Cx, 78.70.Nx, 75.10.Kt, 75.30.Ds 


Geometrically frustrated magnets have been actively inves¬ 
tigated in condensed matter physics illl. In particular, spin ice 
(SI), e.g. (R - Dy or Ho) ||3l3l, provides prototypi¬ 

cal frustrated Ising magnets with the pyrochlore lattice struc¬ 
ture 101, consisting of a three-dimensional network of corner¬ 
sharing tetrahedra [Fig.lHb)]. It displays fascinating features 
such as a finite zero-point entropy and thermally excited 
emergent magnetic or SI monopoles 10,13]. An intriguing the¬ 
oretical proposal for a U(l) quantum spin liquid (QSL) state 
13] has been made for variants of SI endowed with quantum 
spin fluctuations liS. The U(l) QSL state |I3-[I^ is char¬ 
acterized by an emergent U(l) gauge held producing gapless 


fictitious photons and by gapped bosonic spinon excitations 
carrying the SI magnetic monopole charge iiEQ. By 


increasing the transverse interaction, the system can undergo 
a phase transition from the U(l) QSL to a long range ordered 
(LRO) state of transverse spins or pseudospins representing 
electric-quadrupole moments for non-Kramers ions IMl. 
This state can be described as a Higgs phase CSQ. 


In a quest to QSL states in frustrated magneti c sy stems 
from both theoretical 02 ll423n and experimental 0241251] view¬ 
points, an Ising-like pyrochlore Tb 2 Ti 207 (TTO) is a poten¬ 
tial candidate for a U(l) QSL: it has been reported to remain 
in a fluctuating spin state down to 50 mK without magnetic 
LRO 1 2^ 23]. However, the origin of this spin liquid state of 
TTO has been elusive for more than a decade despite many in¬ 
vestigations (see Refs, li, O, 28] and references therein, and 
recent Refs. 123-31]), and is still under hot debate Us, 281. 
To solve this challenging problem of TTO, we start this inves¬ 
tigation by postulating that the theoretically-proposed inter¬ 
action between electric quadrupole moments of non-Kramers 
ions including Tb®^ [the fourth term of Eq. ([D]][il is at work 
for giving the quantum fluctuations to TTO. This postulation 
is a natural consequence of the previous unsuccessful trial- 
and-errors of explaining TTO by taking into account only the 



FIG. 1. (Color online) (a) Phase diagram of the effective Hamilto¬ 
nian Eq. Q determined from CMC simulations. The shading (color) 
represents Tc- Two quadrupole LRO phases, the planar antiferropseu- 
dospin (PAF) and planar ferropseudospin (PF) phases, exist in the 
vicinity of the SI phase JTlil. Classical SI is replaced by a U(l) 
QSL in quantum theory The region enclosed by the dotted line 
represents an acceptable parameter region for the experimental data 
on Tb 2 ,oo 5 Tii, 99507 +j,. The cross mark indicates the typical values 
(S,q) = (0,0.85). (b) Schematic view of the deformation of the /- 
electron charge density due to the PAF order on the pyrochlore lat¬ 
tice. 


interactions between magnetic dipole moments [the first three 
terms of Eq. ([T]i] and the perturbation through first excited 
crystal-field (CE) states 114, and b y taking another as¬ 
sumption of Jahn-Teller (JT) distortion ]^ 33]. Under the 
present postulation, two ground states of off-stoichiometric 
"^ 2 +x^h-xOi+y samples 13411 will possibly be accounted for 
by the U(l) QSL (x < Xc) and electric quadrupolar (x > Xc) 


states of Ref. 

In this Letter, we investigate the hidden order of 
Tb 2 +,tTi 2 -A 07 + 3 , (x = 0.005 > Xc), because the electric 
quadrupolar order is more tractable than the U(l) QSL by us¬ 
ing semi-classical theoretical analyses. Specific heat, magne- 
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tization, and neutron scattering experiments were performed, 
and these experimental data were analyzed using quantum and 
classical Monte Carlo (QMC, CMC) simulations, and a mean- 
field random-phase approximation (MF-RPA). The results 
demonstrate that the hidden order is an electric quadrupolar 
order [Fig.lTfb)] and that the parameters of the model Hamil¬ 
tonian are located close to a phase boundary between the elec¬ 
tric quadrupole and U(l) QSL states [Fig. [TJa)], which sug¬ 
gests that the elusive spin-liquid state of TTOis the U(l) QSL. 
We emphasize that a high-quality single-crystalline sample 
with a well-controlled x value 13511 enables us to accomplish 
this work. 

An effective pseudospin-1/2 Hamiltonian Eh relevant for 
the non-Kramers magnetic doublets of TTO is 


=J„„ ^ cr;cF, ^ e^crj. 

{r.r') r 

e^-e;, 3[e^ • (r - r')][e^, • (r - r')] 


+DrlYj 


r.'|3 


r»/15 


crj.cr;, 


|r-r'|-' |r - r 

+Jnn ^ [2d(o-+o-;, + o-~a-^,) + 0 -^ 0 -^, -i-H.c.)]. 

(T’.T*') 

( 1 ) 


Here, we consider only the CF ground state doublet 3, 
and neglect the first excited doublet at £ ^ 18 K, since we 
are mainly interested in the low-T properties below 2 K. In 
Eq. O, (Tr- are the Pauli matrices (pseudospin) at a site r, 
cr^ = (cr^ + icrl.)l2, and = ±y,0 lElllSl- The mag¬ 
netic dipole moment fisscr^. is parallel to the local (111) axis 
e-; il. The first three terms of Eq. ([T]) represent the nearest- 
neighbor (NN) exchange interaction, the Zeeman energy un¬ 
der a magnetic field H, and the dipolar interaction, respec¬ 
tively. They constitute the classical dipolar SI Hamiltonian 
'Urn 1381] . It can be approximated 13811391] by the NN classi¬ 


cal SI model 'Km.efF = -/nn.eft 2<r,r'> CTrK' “^eft-ff ' Zr 
where 7nn,efr = Jnn + f^nn (f^nn = ^D). The last term of Eq. O 
represents the quadrupole interaction "Kq. We note that the 
transverse components (cr^,ov) of the pseudo spin rep resent 


electric quadrupole (and 16-, 64-pole) moments 111 lLl37l] . Note 


that QMC simulations of the model Eq. ([T]) suffer from a neg¬ 
ative sign problem. On the other hand, thermodynamic prop¬ 
erties away from the QSL state, including phase transitions to 
the LRO phases, can be captured by CMC simulations semi- 
quantitatively. Therefore, in most of the cases, we employ 
CMC simulations. 

CMC simulations were performed up to 1024 pseudospins, 
in which the pseudospin cTr is treated as a classical unit vec¬ 
tor ll4fil] . The resulting zero-field phase diagram is shown in 
Eig.HIa) for the case of Dnn = 0.48 K, namely, fi — 4.6 ijb iH, 
and 7nn = 1.0 K (this value of will be used throughout 
the paper and the choice of 7nn will be explained further be¬ 
low). A quantum mechanical treatment using gauge mean- 
field (ME) theory shows that the classical SI phase region 
in Eig. [TJa) is mostly replaced by a U(l) QSL phase except 
at (5 = q - 0 ill]. The phase diagram has two quadrupole 
LRO phases originating from different ordering patterns of 



FIG. 2. (Color online) (a) Temperature dependence of magnetic 
susceptibility for an applied field along [111]. The solid line is^corr 
(see text for details), (b) Comparison between [.^'exp “TcorrJ * — 
and QMC calculation. Note, 1 emu = 10“^ A m^. 


(cr^,ov): PAE (planar antiferropseudospin) and PE (planar 
ferropseudospin) states denoted in the classical ME phase- 
diagram (Eig. 7 in Ref. d). In particular, a deformation of 
the /-electron charge density iS for the PAE phase is illus¬ 
trated in Eig.lHb) In the following, we will show that 

most of the experimental data on the high-quality polycrys¬ 
talline and single-crystalline samples of TTO with x = 0.005 
can be explained by choosing 7nn = 1 K, 5 = 0, and q - 0.85 
within semi-quantitative analyses. 


Polycrystalline and single-crystalline Tb^+xTi^-i-Oy+v sam- 
ples were prepared by a standard solid-state reaction 13411 and 
by a floating zone method 1^ . Specific heat was measured 
by a quasi-adiabatic method down to 0.1 K using a plate-shape 
crystal with a size of 0.7 x 0.9 x 0.1 mm^ whose shortest di¬ 
mension is along a [110] axis. Magnetization was measured 
by a capacitive Earaday magnetometer using the same sam¬ 
ple. Neutron scattering experiments were performed on NIST- 
BT7 iS using a crystal sample cut from a neighboring part 
of the above sample and on ILL-IN5 i47ll using the powder 
sample with x = 0.005 |34]. 


We first determine the magnitude of 7nn from the magnetic 
susceptibility exp- The T dependence of^exp measured along 
the [111] direction on the single crystal shows an anomaly at 
Tc — 0.53 K [Eig.|2[a)]. While ;^^exp is dominated by the con¬ 
tribution from the CE ground-state doublets, a small but non- 
negligible correction Xcom may arise from higher-energy CE 
states. Thus, we calculated by taking the CE parameters 
of Ref. S and using a single-site approximation, i.e., using 
Eq.(2.L18) of Ref. |^] where the contributions from the CE 
ground state doublet are excluded. The Xcom is also shown 
in Eig. Ea). Now xcndi^ Texp - Tcoit) can be directly com¬ 
pared with a theoretical calculation based on the pseudospin- 
1/2 model Eq. (1). We have performed extensive QMC simu¬ 
lations 1 I 49 I dOT of the nearest-neighbor effective Hamiltonian 
+'?Tq on finite-size clusters up to 1024 pseudospins with 
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FIG. 3. (Color online) (a) Neutron inelastic magnetic spectra of the polycrystalline Tb 2 .005 Til , 99507 +y sample taken at T = 0.1 K. (b) Calculated 
S(Q,E) for the PAF phase with ynn.eif > 0 using Jnn.eir = 1.48 K (y„n = 1.00 K) and {S,q) = (0,0.85). inn,eft is determined by the analysis of 
Xond- (c) Calculated SiQ,E) for the PAF phase with ynn.eif < 0 using Jnn.cir = -1.77 K (y„n = -2.25 K) and (6,q) - (-0.5, -1.0). 7nn,eff is 
determined by the analysis of^ond- (d) T -dependence of intensities of the single-crystal neutron Bragg scattering at (002) and ( 5 1|). 


typical Monte-Carlo steps of 200000. The experimental data 
OfGnd) are well reproduced by the QMC calculations in a wide 
range of 6 and q, if we take |Tnn,effl = 1.3 - 1.9 K. Note that 
because of the negative sign problem of the QMC simulation, 
the analyses have been limited to a relatively high tempera¬ 
ture range, 5 <T < 15 K. In Fig.|3b), we show a representa¬ 
tive comparison between xcnd and the QMC results obtained 
for (b, q) - (0,0.85) as determined below. This comparison 
yields Tnn.eff = 1.48(1), leading to 7nn = 1.0(1) K. This value 
of Tnn.eff Is of the Same order as the previous estimation 115111 . 

Next, we confirm the positive sign of Tnn.eft and extract the 
parameter values of (b, q) from the comparison of the inelas¬ 
tic neutron scattering data. The previous inelastic magnetic 
neutron-scattering spectra measured on the x = 0.005 powder 
sample at T = 0.1 K (<s: Tc), have shown a nearly flat broad 
peak at 0.1 meV in the (Q, E) space ll34ll . as shown in Fig. 3(a). 
The peak is broader than the instrumental resolution, which 
suggests dispersive excitations. This experimental behavior 
can be described in terms of pseudospin waves in the PAF and 
PF phases, as discussed in Ref. 13711 . The powder-averaged 
dynamical magnetic structure factor ^(IQI,^) is calculated 
within the MF-RPA Q, which can correctly describe the 
spectrum in an ordered state within the linear spin-wave ap¬ 
proximation. Extensive calculations in a wide range of the pa¬ 
rameters (b, q) show that reasonable agreements are obtained 
in the PAF phase with (b, q) = (0.0 + 0.4,0.8 + 0.3) [Fig.I^b)] 
and in the PF phase with (6,q) - (-0.54 + 0.02,q < 1.1), 
when we fix = 1.0(1) K [Fig. [11^)]. We note that only 
the cases of Tnn.eft > 0 can reasonably reproduce the observed 
features for the case of the PAF phase [Fig. [3b)]. The case of 
Tnn.eft < 0 gives highly dispersive spectra that are not compat¬ 
ible with the experimental results [Fig.|3c)]. 

Finally, we show, using CMC simulations of Eq. (1), that 
the parameter set showing the PAE explain reasonably well 
the observed specific heat Cp{T,H) under weak [111] held 


while the other parameter set showing the PE does not. Eig- 
ure HJa) shows the T dependence of Cp{T,H) under [111] 
held up to 1.0 T. The sharp peak at survives only up to 
fioH - 0.1 T, turning into broad double peaks at 0.3 T. A 
full map of Cp{T,H) is shown in Eig.|4jb). Eor comparison, 
maps of the calculated specific heat C{T, H) by CMC simula¬ 
tions are presented in Eigs. |4jc) and (d) for the same param¬ 
eter choices as determined above, namely, (b, q) - (0,0.85) 
and (-0.54,0.5), respectively. Clearly, the PAE case shows 
a better qualitative agreement with the experiment, although 
with some discrepancy in the magnetic held and temperature 
scales. We note that CMC simulations with these PAE param¬ 
eters also reproduce the experimental results of C{T, H) under 
the [100] field i^. 

One may not clearly see the change of states under the [111] 
held below 0.4 K in Cp(T,H) [Eigs. EJa) and (b)]. How¬ 
ever, certain changes are observed in the magnetization M. 
Eigure [Sja) shows M-H curves at several temperatures un¬ 
der the [111] held. Two clear step-like kinks are observed 
at ij.qH\ ^ 0.14 T and ;Uo //2 - 1.3 T below the zero-field 
Tc. CMC simulations with the same parameters as used 
above demonstrate the first kink at ^qH\ [Eig|3b)], indicat¬ 
ing that it is a crossover or a phase transition from the three 
dimensional (30) PAE state 191] to the two dimensional (2D) 


PAE state 11 ll l40l] . This means that in intermediate fields 
(/io/f - 0.5 T), the system behaves as decoupled 2D kagome 
layers of quadrupole moments separated by triangular layers 
of polarized magnetic moments. This bears resemblance to 
the kagome ice state of SI materials ll53n . In contrast, the sec¬ 
ond kink appears at a higher held than that of experiments at 
IJ.qH 2 . This result suggests that higher order terms neglected 
in Eq. (HJ, such as terms due to CE excited states iItHi . are re¬ 
quired for further explanation of the behavior at fields higher 
than about 1 T. 

All the above comparisons between the experiments and 
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FIG. 4. (Color online) (a) Temperature dependence of the observed specific heat Cp(T,H) for H || [111], (b) Temperature-field map 
of Cp(T,H) for H || [111]. Filled circles in the map are peak positions of Cp(T,H). (c,d) CMC results of specific heat C(T,H) for (c) 
(S,q) = (0,0.85) with = 1.00 K and for (d) {6,q) = (-0.54,0.5) with = 0.92 K. The values of have been determined from the 
comparison of xond with the QMC results for each case. Solid, dashed, and dotted lines in (b), (c), and (d) are guides to the eyes. Labels 
in maps (b) and (c) indicate assigned states from the analysis with Eq.JTll and CMC simulations Hi; i.e., (I) a paramagnetic paraquadrupole 
state, (II) the 3D PAF state, and (III) the 2D PAF state. 



FIG. 5. (Color online) Field dependence of the magnetiza¬ 
tion M under the [111] magnetic field, (a) Experimental data on 
Tb 2 .oo 5 Tii, 99507 + 3 , above and below Tc = 0.53 K. (b) CMC results 
obtained with the same parameters as for Figs.[3tb) and|4lc). 


theories show that the LRO of TTO is the quadrupole or¬ 
der characterizing the PAF phase. Although this LRO can¬ 
not be detected directly through neutron Bragg scattering, we 
observed some indication of it. In fact, weak magnetic re¬ 
flections are observed at the forbidden ( 002 ) position and at 
the superlattice ( 555 ) position [Fig. [3d)]. Polarized neutron 
scattering experiments at BT7 confirm that both of them are 
magnetic. The long-range ordered magnetic moments of these 
reflections are roughly ~ 0.1 fis, which is too small to be the 
primary order parameter. We speculate that the (002) reflec¬ 
tion appears simultaneously with the PAF order, whose order 
parameter is characterized by the wave vector k - 0 11113], 
and is induced by higher order terms neglected in Eq. (IB- On 
the other hand the ( 555 ) reflection, observed also in a powder 
sample iHI] . suggests a different origin because of the differ¬ 
ent T -dependence for the (002) reflection [Fig. [3d)]. 


Electric quadrupolar orders are related to the deformation 
of /-electron charge densit y. T hese naturally couple to dis¬ 
placements of ligand ions OTSll and may induce cooperative 
JT effects and JT structural distortions [33, H]. Since the 


quadrupole-coupling terms of Eq. ([B are derived as the elec¬ 
tronic coupling but are symmetry allowed terms under the 
space group of the pyrochlore lattice El , these may contain 
a phonon-coupling contribution illsl] . Thus, a direct detec¬ 
tion of the quadrupole order (using resonant X-ray scattering) 
and/or of a small associated JT lattice distortion are difficult 
yet interesting topics for future investigations. 

In summary, the hidden order of Tb 2 +ATi 2 -j: 07 +y with x = 
0.005 > Xc has been studied using thermodynamic and neu¬ 
tron scattering measurements on single-crystalline and poly¬ 
crystalline samples under the control of the off-stoichiometry 
of X. We take account of magnetic-dipole and electric- 
quadrupole moments of the CE ground-state doublet of the 
non-Kramers Tb^^ ion as well as the theoretically proposed 
quadrupole interaction El- Semi-quantitative analyses of the 
experimental data based on a simple pseudospin-1/2 Hamilto¬ 
nian demonstrate that the hidden order is an order of the elec¬ 
tric quadrupole moments [Fig. [Bb)]- The estimated model 
parameters are located close to the phase boundary between 
the quadrupolar and U(l) QSL states. This result implies that 
the putative SL state of TTO studied for more than a decade is 
)SL. Investigations in the context of a Higgs transi- 
llOl l56l] and on the relation of it with the neighboring 
U(l) QSL phase using Th 2 +iJ^'i 2 -xOi+y single crystals 0351] are 
fascinating future topics. 
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Supplemental Material: Quadrupole Order in the Frustrated Pyrochlore T\)2+x^i2-xOi+y 

ABSTRACT 


In this supplemental material, we discuss effects of higher excited crystal held (CF) states of Tb 2 Ti 207 (TTO) using the 
analysis of the magnetization M. We also present the result of inelastic neutron scattering (INS) data above and discuss a 
local symmetry reduction which couples to the static Jahn-Teller distortion but is thought to be a very minor effect for the present 
system. 


EFFECTS OF HIGHER EXCITED CF STATES: M-H CURVE UNDER [111] FIELD 


When we consider TTO in the wider T and wider H ranges beyond the scope of the present investigations, we cannot neglect 
the contribution from the hrst excited CF doublet state (E ^ 18 K) & An interesting quantum effect coming from the hrst 
excited CF state was pointed out in Ref. llS2|] . The NN exchange couplings between the magnetic dipole moments (J, ■ Jj) gives 
rise to a perturbative term to the effective S - 112 Hamiltonian in the same form as the 6 0 term in Eq. (1) of the main text. 

Thus the effective pseudospin S - 112 Hamiltonian that we postulate in this work (Eq. (1) in the main text) may be regarded as 
a renormalized 5 = 1/2 Hamiltonian including effects of the higher excited CE states. 

In this section, we try to explain how the higher excited CF states affect the present system from an experimental viewpoint. 
To this end the magnetization curve under the [111] magnetic held is a good example. A magnetization curve within a single-site 
approximation can be calculated using the single-site CF Hamiltonian |S3t]. We thus used the CF parameter set given in Ref. 11541] 
and took account of the four crystallographic sites. 

The M-H curve at 0.1 K is calculated and plotted in Fig.[S2] In this hgure, we also plotted the experimentally observed M-H 
curve at 0.12 K (the same data in Fig. 5(a) of the main text) for comparison. The obvious discrepancy between these two M-H 
curves should be accounted for by multi-site effects, i.e., by quantum or classical many-body effects of this frustrated system. 
Note that the gradual increase of the calculated M in ^qH > 0.5 T originates mainly from the hybridization of the hrst excited 
CF states to the CF ground states. This is a simple single-site quantum effect. A similar gradual increase is also observed in the 
experimental M-H curve. It is thus considered that a useful H range of the effective S - 112 Hamiltonian (Eq. (1) of the main 
text) is < 1 T. 



FIG. SI. (Color online) Magnetization curves under the [111] magnetic held at 0.1 K. The calculated M-H curve within the single-site 
approximation and the experimentally observed M-H curve at 0.12 K are shown. 


In Eig. ISll the calculated M-H curve exhibits a quick rise to M = 2.6 /ib in EqH < 0 05 T. This is naturally considered as that 
the magnetic dipole moment of Tb^'^ is fully polarized along a local [111] axis within the CE ground state doublet in low held. 
On the other hand, the experimental M-H curve exhibits a slower increase around ^qH = 0.1 T and the kink at ^qH\ = 0.14 T 
(M =1.2 yUfi/Tb). The classical Monte Carlo (CMC) simulation shown in Eig. 5(b) of the main text approximately reproduces 
these experimental behaviors. We can thus think a classical many-body effect is at work around the kink to certain extent. 
Around the hrst kink of fi^Hi for the CMC simulation, the magnetization is about M = 1.5 yUfi/Tb = (4.6/ 3) ;UB/T b. This 
implies that the state bears certain resemblance to the 1/3-plateau state or the kagome ice state of SI materials iS5llS6l] . In this 
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state, while the magnetic moments in the triangular layers are polarized, the electric quadrupole moments in the kagome layers 
behave as a 2D system. It can be regarded as a 2D frustrated system of quadrupole moments. 

Here we also note that CMC simulations qualitatively reproduce the temperature dependence of the magnetic susceptibility 
along the [111] direction 10 . It shows a similar anomaly at to that observed in experiments (Fig. 2(a) of the main text). This 
result supports the conclusion that Tb 2 +;cTi 2 -A 07 + 3 , can be basically represented by the simple effective S = 1/2 Hamiltonian 
(Eq.(l) of the main text), although there could be non-negligible and complicated terms for a full description. 


TREATMENT OF THE MAGNETIC MOMENT OF THE GROUND STATE DOUBLET 


In this section, we describe the treatment of the magnetic moment of the ground state doublet ^ through the calculations and 
its effect to the estimation of parameters in Eq. (1) of the main text. It is known that there are several reports on the esti matio n 
of for TTO: = 5.1 (Ref. 0 ), 4.8 (Ref. 0 ), 3.9 (Ref. 0 ), and 5.4 fiB/T'h^'^ (Ref. lHH). 

Therefore, there is some ambiguity in the choice of //. This fact also affects the choice of Dnn for the calculations. In this context, 
we adopted yu = 4.6 in the calculations since it reproduces reasonably well the observed value of M at the first kink of 

jUoHi. However, even when using other choices of yu, the main conclusion of the present study does not alter. It only changes 
parameter values of J^n, 6, and q slightly; for example, if we use ju = 5.3 yUB/Tb^^ (i.e.. Dm, = 0.63 K), we obtain 7nn = 0.9(1) K, 
and (6,q) - (0.0 + 0.4,1.0 + 0.3) for reproducing experimental results. With these parameter values, we still obtain the PAE 
phase, although the boundary between the SI and PAE phases should be changed slightly from that shown in Eig. 1(a) of the 
main text. 


INELASTIC NEUTRON SCATTERING SPECTRUM AT TEMPERATURES ABOVE 


In this section, we present the INS data of the Tb 2 . 005 Tii. 995 O 7 + 3 , powder sample above and discuss a possibility of a local 
symmetry reduction which affects the splitting of the INS spectrum. 

Eigure|S2] shows the INS spectrum at T = 0.7 K. A quasielastic scattering spectrum is observed. This spec trum changes into 
an inelastic scattering spectrum with decreasing temperatures below Tc (see Fig. 3(a) of the main text and Ref. iH). This result 
indicates that the local symmetry reduction induced by the substitution of an element like x of Tb 2 +ATi 2 _A 07 + 3 , is a very minor 
effect or is absent for the case of Tb 2 . 005 Tii. 995 O 7 +v. This is because the local symmetry reduction often inevitably couples to 
the static Jahn-Teller effect which may lead to the splitting of the INS spectrum persisting to temperatures above T^. Therefore, 
the result of Fig. [S2] indicates that the splitting of the peak (Fig. 3(a) of the main text) is mainly induced by another effect. We 
discussed it by considering the development of a pseudospin wave in the quadrupolar LRO state and demonstrated a similar 
spectrum by MF-RPA calculations (see the main text and Ref. ||S12|] for the details). 
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FIG. S2. (Color online) Neutron inelastic magnetic spectrum at T = 0.7 K (above D = 0.53 K). 
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